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Abstract 

The detection of the theoretically expected dark matter is central to particle 
physics and cosmology. Current fashionable supersymmetric models provide a natu- 
ral dark matter candidate which is the lightest supersymmetric particle (LSP). Such 
models combined with fairly well understood physics like the quark substructure of 
the nucleon and the nuclear structure (form factor and/or spin response function), 
permit the evaluation of the event rate for LSP-nucleus elastic scattering. The thus 
obtained event rates are, however, very low or even undetectable. So it is imperative 
to exploit the modulation effect, i.e. the dependence of the event rate on the earth's 
annual motion. Also it is useful to consider the directional rate, i.e its dependence 
on the direction of the recoiling nucleus. In this paper we study such a modulation 
effect both in non directional and directional experiments. We calculate both the dif- 
ferential and the total rates using both isothermal, symmetric as well as only axially 
asymmetric, and non isothermal, due to caustic rings, velocity distributions. We 
consider We find that in the symmetric case the modulation amplitude is small. The 
same is true for the case of caustic rings. The inclusion of asymmetry, with a real- 
istic enhanced velocity dispersion in the galactocentric direction, yields an enhanced 
modulation effect, especially in directional experiments. 



1 Introduction 

In recent years the consideration of exotic dark matter has become necessary in 
order to close the Universe Q' ||. Furthermore in in order to understand the large 
scale structure of the universe it has become necessary to consider matter made 
up of particles which were non-relativistic at the time of freeze out. This is the 
cold dark matter component (CDM). The COBE data @ suggest that CDM is at 
least 60% Q. On the other hand during the last few years evidence has appeared 
from two different teams, the High-z Supernova Search Team || and the Supernova 
Cosmology Project > which suggests that the Universe may be dominated 
by the cosmological constant Ar. As a matter of fact recent data the situation 
can be adequately described by a barionic component £Ib = 0.1 along with the 
exotic components Qcdm = 0.3 and f?A = 0.6. In another analysis Turner []|] gives 
Q m = Ucdm + m B = 0.4. Since the non exotic component cannot exceed 40% of 
the CDM ||' Q, there is room for the exotic WIMP's (Weakly Interacting Massive 
Particles). In fact the DAMA experiment [|l0| has claimed the observation of one 
signal in direct detection of a WIMP, which with better statistics has subsequently 



been interpreted as a modulation signal [|TT | . 



The above developments are in line with particle physics considerations. Thus, 
in the currently favored supersymmetric (SUSY) extensions of the standard model, 
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the most natural WIMP candidate is the LSP, i.e. the lightest supersymmetric 
particle. In the most favored scenarios the LSP can be simply described as a Majo- 
rana fermion, a linear combination of the neutral components of the gauginos and 
Higgsinos |]> 

Since this particle is expected to be very massive, m x > 30GeV , and extremely 
non relativistic with average kinetic energy T < WOKeV , it can be directly detected 



[12] [13] mainly via the recoiling of a nucleus (A,Z) in the elastic scattering process: 

x + (A,Z) - x + (A,Z)* (1) 

(x denotes the LSP). In order to compute the event rate one needs the following 
ingredients: 

1) An effective Lagrangian at the elementary particle (quark) level obtained in 
the framework of supersymmetry as described in Refs. 0], Bottino et al. pO] and 



2) A procedure in going from the quark to the nucleon level, i.e. a quark model 
for the nucleon. The results depend crucially on the content of the nucleon in quarks 
other than u and d. This is particularly true for the scalar couplings as well as the 
isoscalar axial coupling 111, EH] . 



3) Compute the relevant nuclear matrix elements [|2g]~[[30|] using as reliable as 
possible many body nuclear wave functions. By putting as accurate nuclear physics 
input as possible, one will be able to constrain the SUSY parameters as much as 
possible. The situation is a bit simpler in the case of the scalar coupling, in which 
case one only needs the nuclear form factor. 

Since the obtained rates are very low, one would like to be able to exploit the 
modulation of the event rates due to the earth's revolution around the sun. To 



this end one adopts a folding procedure assuming some distribution [g, ^T], 32] of 
velocities for the LSP. One also would like to know the directional rates, by observing 
the nucleus in a certain direction, which correlate with the motion of the sun around 
the center of the galaxy. 

The purpose of our present review is to focus on this last point along the lines 
suggested by our recent work [JD], [l8|]. For the reader's convenience, however, we 
will give a description in sects 2, and 4 of the basic SUSY ingredients needed to 
calculate LSP-nucleus scattering cross section. We will not, however, elaborate on 
how one gets the needed parameters from supersymmetry. The calculation of these 
parameters has become pretty standard. One starts with representative input in 
the restricted SUSY parameter space as described in the literature, e.g. Bottino et 



al. [20], Kane et al., Castano et al. and Arnowitt et al. |nj[. Our own SUSY input 



parameters will appear elsewhere [19] 

After this we will specialize our study in the case of the nucleus 127 1, which is 
one of the most popular targets []l(|-|33| |34|]~ [36]. To this end we will consider 
velocity distributions both isothermal, symmetric Maxwell-Boltzmann [g] as well as 
only axially asymmetric, like that of Drukier |3^] and non isothermal as well. In the 
non isothermal case we will consider the relatively simple case in which dark matter 
accumulated as late in-fall of such matter into our galaxy, i.e. the Sikivie caustic 
rings fT|. 

Since the expected rates are extremely low or even undetectable with present 
techniques, one would like to exploit the characteristic signatures provided by the 
reaction. Such are: a) The modulation effect, i.e the dependence of the event rate 
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on the velocity of the Earth and b) The directional event rate, which depends on the 
velocity of the sun around the galaxy as well as the the velocity of the Earth. The 
latter effect, recognized sometime ago |37| has recently begun to appear feasible by 
the planned UKDMC experiment We will study both of these effects in the 

present work. 

In all calculations we will, of course, include an appropriate nuclear form factor 
and take into account the influence on the rates of the detector energy cut off. We 
will present our results a function of the LSP mass, m x , in a way which can be 
easily understood by the experimentalists. 



2 The Nature of the LSP 



Before proceeding with the construction of the effective Lagrangian we will briefly 
discuss the nature of the lightest supersymmetric particle (LSP) focusing on those 
ingredients which are of interest to dark matter. 

In currently favorable supergravity models the LSP is a linear combination j|, 12] 
of the neutral four fermions B, W3, H\ and H2 which are the supersymmetric part- 
ners of the gauge bosons and and the Higgs scalars Hi and H2. Admixtures 
of s-neutrinos are expected to be negligible. 

In the above basis the mass-matrix takes the form pi [23| 

/ 



V 



Mi 


-m z Cj3S w 
m z Sj3S w 





M 2 

m z Cj3C w 
-m z C/3C u 



-m z cps w 
m z Cj3C w 


-A* 



m z sps w \ 
-m z S[ic w 
-A* 
/ 



(2) 



In the above expressions cw = cosOw, sw = sinOw, eg = cos/3, sp = sin/3, 
where tan/3 = {v2)/{v\) is the ratio of the vacuum expectation values of the Higgs 
scalars H2 and H\. u is a dimensionful coupling constant which is not specified by 
the theory (not even its sign). The parameters tan/3, M\, M2, u are determined by 



the procedure of Kane it et al and Castano et al in Ref. [21| using universal masses 
of the GUT scale. 

By diagonalizing the above matrix we obtain a set of eigenvalues mj and the 
diagonalizing matrix Cij as follows 



( Br \ 




( XlR \ 




( Bl \ 




1 XlL \ 


W 3 R 
HlR 


= (eg) 


X2R 
X3R 








X2L 
X3L 


V H 2R j 




\ X4R J 




\ H2L ) 




V X4L J 



(3) 



with C% 



CI -e lX i The phases are Aj 



, . ,„... „, ,,, - 0, 7T depending on the sign of the eigemass. 

Another possibility to express the above results in photino-zino basis 7, Z via 



W3 = sinOwl — cosOwZ 
Bq = cosOwl + sin6y/Z 



(4) 



In the absence of supersymmetry breaking {M\ = M2 = M and u = 0) the photino 
is one of the eigenstates with mass M. One of the remaining eigenstates has a 
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Table 1: The essential parameters describing the LSP and Higgs. For the definitions see 
the text. 



Solution 


#1 


#2 


#3 


#4 


#5 


#6 


#7 


#8 


#9 


m x (GeV) 


126 


27 


102 


80 


124 


58 


34 


35 


50 


m h 


116.0 


110.2 


113.2 


124.0 


121.0 


105.0 


103.0 


92.0 


111.0 


m H 


345.6 


327.0 


326.6 


595.0 


567.0 


501.0 


184.0 


228.0 


234.0 


m A 


345.0 


305.0 


324.0 


594.0 


563.0 


497.0 


179.0 


207.0 


230.0 


tan2o 


0.245 


6.265 


0.525 


0.410 


0.929 


0.935 


0.843 


1.549 


0.612 


tan/3 


10.0 


1.5 


5.0 


5.4 


2.7 


2.7 


5.2 


2.6 


5.3 



zero eigenvalue and is a linear combination of H\ and H 2 with mixing sin/3. In 
the presence of SUSY breaking terms the B, W3 basis is superior since the lowest 
eigenstate xi or LSP is primarily B. Prom our point of view the most important 
parameters are the mass m x of LSP and the mixings Cj\,j = 1,2,3,4 which yield 
the xi content of the initial basis states. These parameters which are relevant here 
are shown in Table 1. We are now in a position to find the interaction of xi with 
matter. We distinguish three possibilities involving Z-exchange, s-quark exchange 
and Higgs exchange. 



3 The Feynman Diagrams Entering the Direct 
Detection of LSP. 

The diagrams involve Z-exchange, s-quark exchange and Higgs exchange. 

3.1 The Z-exchange contribution. 

This can arise from the interaction of Higgsinos with Z which can be read from Eq. 



C86 of Ref. [23 1 



L = -^-\[H lRltl H 1R - HilI^Hil ~ (H 2R ^H 2R - H 2L ^H 2L )]Z» (5) 

COSUW 4 

Using Eq. (^) and the fact that for Major ana particles xi^X = 0, we obtain 

L = ^r-ld^il 2 - |C 4 i| 2 )xi7m75Xi^ (6) 
cost>w 4 

which leads to the effective 4-fermion interaction 
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where the extra factor of 2 comes from the Majorana nature of Xi- The neutral 



hadronic current J? is given by 



J\ = -<H\{-^sin 2 e w - ~(1 - 75) - sin 2 6 w r 3 }q (8) 
at the nucleon level it can be written as 

Jf = -iV7 A { sin 2 6 w - gv(^~ sin 2 6 w )T 3 + ^gA75n}N (9) 



Thus we can write 



L eff = -9z {xa ^ Xl)Jx{Z ) (10) 



where 
and 



J X (Z) = N lx [f v {Z) + f v {Z)r 3 + fl(Z) l5 + /A(Z) 75 r 3 ]iV (11) 



m 



2 



f v {Z) = 2(\C 31 \ 2 -\C 41 \ 2 ) ZI f±_ sin 2 e w (12) 



m 



ti{Z) = -2(|C 3 i| 2 -|C4i| 2 ) ^ 2 9v(^-sin 2 9 w ) (13) 

m z — q z 2 

f°A(Z) = (14) 

f\(Z) = 2 {\C 3l \ 2 -\C Al \ 2 )-^- 2 \g A (15) 

with gy = 1.0, gA = 1-24. We can easily see that 

WK = -»^-i)=-ui m 

Note that the suppression of this Z-exchange interaction compared to the ordinary 
neutral current interactions arises from the smallness of the mixings C31 and C41, 
a consequence of the fact that the Higgsinos are normally quite a bit heavier than 
the gauginos. Furthermore, the two Higgsinos tend to cancel each other. 



3.2 The s-quark Mediated Interaction 

The other interesting possibility arises from the other two components of xii namely 
B and W 3 . Their corresponding couplings to s-quarks can be read from the appendix 
C4 of Ref. |23j They are 

L eff = -gV2{q L [T 3 W m -tane w (T 3 -Q)B R ]q L 

- tan9 w q R QB L q R } + HC (17) 

where q are the scalar quarks (SUSY partners of quarks). A summation over all 
quark flavors is understood. Using Eq. (||) we can write the above equation in the 
Xi basis. Of interest to us here is the part 

L eff = gV2{(tane w (T 3 -Q)C^-T 3 Ci)q LX iRqL 

+ tan9 w C n Qq R xiLQR} (18) 
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The above interaction is almost diagonal in the quark flavor. There exists, however, 
mixing between the s-quarks qi and q R (of the same flavor) i.e. 

qi = cosOqqi + sinOqfe (19) 

q R = -sinOqqi + cosOgfc (20) 

with 



m 



m u R + m 2 cos2(3/2 



d - m 2 - + m 2 z cos2P/2 

Thus Eq. becomes 

I e // = gV2{[B L cos6q q L xiR~ B R sin9qq R xiL}qi 
+ [B L sin8q-q L xiR + B R cos9qq R XiL] fa] 

with 

Bl(q) = —\Ciitan9 u — ^C^, g = u {charge 2/3) 

Bl(o) = -\C$tan0u + |C^, g = d {charge - 1/3) 

Br(<1) = %tanO u Cu, q = u {charge 2/3) 

B R {q) = — jjianfl^Cii, g = <i {charge — 1/3) 

The effective four fermion interaction takes the form 

L e// = {gV2) 2 {{B L cos9qq L xi R - B R sin9qq R xiL) 
~2— {B L cos9 q xi R qL - B R sin9qXiLqn) 

q m qf 

1 

+ {B L sin9 q q L xi R + cos9q-q R XiL) 



{B L sin9 q xi R qL + B R cos9qXi L q R )} (23) 

The above effective interaction can be written as 

Leff = L L e f+ RR + L L e R (24) 

The first term involves quarks of the same chirality and is not much effected by 
the mixing (provided that it is small) . The second term involves quarks of opposite 
chirality and is proportional to the s-quark mixing. 

i) The part L^ f +RR 
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Employing a Fierz transformation L^j~ RR can be cast in the more convenient 
form 



L%+ RR = ( 9 V2?2{-\){\B L ? 



cos 2 9g sin 2 bq _ A 



i 2 — rriz2 



+ —mixqLXlRl XiR 



?1 ^ ""95 



sin 



2a_ „„„2z 



+\B R \\- q - + q -)q R lxqRXiLl\iL} (25) 

q z — ma. q z — m-,2 

The factor of 2 comes from the Majorana nature of LSP and the (-1/2) comes from 
the Fierz transformation. Equation (^) can be written more compactly as 

L e ff = -—^2{q7\(PoR + (3mT3)(l + 7s)<? 



qi\{PoL + P 3 lt 3 )(1 - 75M(X17 7 Xi} (26) 



with 



Por = (^xl R + ^X 2 d R )\Cntan6 w \' 



PsR = {lxl R -lx} R )\Cntane w \ 2 (27) 
Pol = \lc R tane w +lc R \ 2 xl r +\lc R tane w -lc R \\} 



d L 



hi = \lc R tan0 w + lc R \ 2 xl L -\lc R tan0 w -^C R \ 2 X } L 



with 



2 2 
2 2 m W , 2 m VK 

XgL ~ C «m d 2 -q 2 S «m d 2 - q 2 
i\ 12 

2 2 
2 2 m VK 1 2 m W 

XqR °5 9 ~r ^5 o 

y * — q z H m-,2 — q z 

1l ^2 
Cq = COsOg, Sq = sinOq (28) 

The above parameters are functions of the four-momentum transfer which in our 
case is negligible. Proceeding as in Sec. 2.2.1 we can obtain the effective Lagrangian 
at the nucleon level as 

L™ = -^(xi7 A 7 5 Xi)Ja(^) (29) 



with 



Jx{q) = N lx {ft{q) + ti(q)n + f A (q)j 5 + fXiqH^N (30) 

Sv = 6 (A)ii - Pol), ft = 2(Am - Pzl) 
f°A = 2g° a g v (PoR + Pol), A = ^9a(P S r + Psl) (31) 

with g v = 1.0 and g A = 1.25. The quantity g A depends on the quark model for the 
nucleon. It can be anywhere between 0.12 and 1.00. 
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We should note that this interaction is more suppressed than the ordinary weak 
interaction by the fact that the masses of the s-quarks are usually larger than that of 
the gauge boson Z°. In the limit in which the LSP is a pure bino {C\\ = 1, C21 = 0) 
we obtain 



2 / 4 
(3 m = tan Qw[g 

tan 2 w , 2 , 2 n 
Pol = 36 UtL+XiJ 



4 2 1 2 

^~ U R g ^d R 



tan 2 9 w , 2 2 \ / QO \ 
&l = 36 (Xfi^-XjJ (32) 

Assuming further that Xu R = Xi R = Xu L = Xg L we obtain 

* +^ (33) 

If, on the other hand, the LSP were the photino (Cn = cosOw-, C21 = sinOw, C31 = 
C41 = 0) and the s-quarks were degenerate there would be no coherent contribution 



(/£ = if /3 L = Am) 



From Eq. (23) we obtain 

L L p f = -{g^) 2 sin29~ q B L {q)B R {q)\' ' 



ff — ~\yv * o^'^g-^LKI ^>H\1 7; ~2 2 

3 J 2 q z — rnsfl q A — m~2 

ii 12 

(qixiRXiLqa + mxiLXinqi) 

Employing a Fierz transformation we can cast it in the form 

G F 

L e ff = -—^[P+(q)(qqxixi + 975^x175x1 - (g<v?)(xio" |W 'xi)) 

+P-{qnqx\Xi + <F375<?Xi75Xi ~ ^V^Xi^Xi)] 
where for the light quarks u and d 

P± = -tan6 w C n {2sin29u{-C$tan9 w + -C^]Aq 
3 b 2 

[-CfitanOw - -< 
b 2 

for quarks other than u and d we only have only the isoscalar contribution which is 
given by 

2 11 

f3 + = -tan6 w Cu{2sin29u[-CRutan9w + -CR 2 i]A il 

3 b 2 



T sin29£c^tane w - -C&]Aj} 



+ sin29 3 [^CR u tan6 w - ^CfjAj} 



[-CRntan0 w - ]-( 
b 2 

Where in the last expression u indicates quarks with charge 2/3 and d quarks with 
charge -1/3. In all cases 
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and an analogous equation for Aj. 

The appearance of scalar terms in s-quark exchange has been first noticed by 
Griest. [16| It has also been noticed there that one should consider explicitly the 



effects of quarks other than u and d [14] in going from the quark to the nucleon 



level. We first notice that with the exception of t s-quark the q~L — Qr mixing small. 
Thus 



sin29,r,Au 



sin29jAd 



2m u (A + ucot^m^r 
(m\ L - q 2 ){m\ R - q 2 ) 

2m ( i(A + ntan^m^ 



d (m 2 - — « 2 )(m 2 T — a 2 ) 

In going to the nucleon level and ignoring the negligible pseudoscalar and tensor 
components we only need modify the above expressions for all quarks other than t 
by the substitution m q — > / g mjv- We will see in the next section that the quarks s,c 
and b tend to dominate. For the t s-quark the mixing is complete, which implies 
that the amplitude is independent of the top quark mass. Hence in the case of 
the top quark we do not get an extra enhancement in going from the quark to the 
nucleon level. In any case this way we get 



L eff = ^[f°(q)NX +/i(7)AY,A']viVi (•">.!! 

with 



= /?/?+ and fl(q) = flf3- (35) 

(see sect. 3.3 for details). In the allowed SUSY parameter space considered in this 
work this contribution can be neglected in front of the Higgs exchange contribution. 
This happens because for quarks other than t the s-quark mixing is small. For the 
t-quark, as it has already been mentioned, we have large mixing, but we do not get 
the advantage of the mass enhancement. 



3.3 The Intermediate Higgs Contribution 

The coherent scattering can be mediated via the intermediate Higgs particles which 
survive as physical particles . The relevant interaction can arise out of the Higgs- 
Higgsino-gaugino interaction which takes the form 

l h X x = ^ {WrH^lHI * - W R H XL Hl * 

- tane w {B R H 2L H$*-B R H 1L H?*j)+H.C. (36) 



Proceeding as above we can express W an B in terms of the appropriate eigenstates 
and retain the LSP to obtain 

L = -^{{C^-tan6 w C^)C A1 xiRXiLH 2 * 

- (Ci-tan9 w C^)C 31 xiRXiLHr)+H.C. (37) 
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We can now proceed further and express the fields H®* , H®* in terms of the 
physical fields h, H and A. The term which contains A will be neglected, since it 
yields only a pseudoscalar coupling which does not lead to coherence. 

Thus we can write 

C eff = -^xxN[fs(H) + fg(H)T 3 )N (38) 

where 

fs(H) = -(g u + 9d)+9s + gc + g b + gt (39) 

fl(H)= l -{g u -g d ) (40) 

with 

9a t = [gi(h)— -3+52 (H — - , a i = n,c,t 41 
smp smp rriN 

,,.sina , rr .cosa,m K . , , , . 

5 Ki = l-9l(h) -+ 92(H) -]-*, K t = d,s,b 42 

cosp cosp m-N 

gi (h) = 4(CEtan9 w - C^)(C 41 cosa + C 31 sina) m " mw ^ 

m h — q l 

g 2 (H) = A(C?vbanO w - CSVC^sina - C 31 cosa) m ^ mw (44) 

m H ~ 1 

where mjy is the nucleon mass, and the parameters nih, ran and a depend on the 
SUSY parameter space (see Table 1). If one ignores quarks other than u and d 
(model A) and uses m u = 5MeV = m d j2 finds [24] 

/ s ° = 1.86( 5u + 9d)/2, fl = 0.49( 5m - g d )/2, (45) 



4 Going from the Quark to the Nucleon Level 



As we have already mentioned, one has to be a bit more careful in handling quarks 



other than u and d since their couplings are proportional to their mass [14]. One 
encounters in the nucleon not only sea quarks (uu, dd and ss) but the heavier quarks 
also due to QCD effects [15). This way one obtains the scalar Higgs-nucleon coupling 



by using effective quark masses as follows 

f u m N , m d -> fdm N . 



f s m N 



yrriN, (heavy quarks c,b,t) 



where rriN is the nucleon mass. The isovector contribution is now negligible. The 
parameters f q , q = u,d,s can be obtained by chirala symmetry breaking terms in 
relation to phase shift and dispersion analysis. Following Cheng and Cheng [^] we 
obtain 



Su = 0.021, Sd = 0.037, Ss = 0.140 



(model B) 
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/„ = 0.023, f d = 0.034, f s = 0.400 (model C) 

We see that in both models the s-quark is dominant. Then to leading order via 
quark loops and gluon exchange with the nucleon one finds: 

f Q = 2/27(1- £,/,) 

This yields: 

f Q = 0.060 (model C) 

f Q = 0.040 (model C) 

There is a correction to the above parameters coming from loops involving s-quarks 
[15]. The leading contribution can be absorbed into the definition if the functions 
gi(h) and g2(H) as follows : 

sin(a+0) sinf}\ 



,2 

_l : — 

cos z 0w s%na> 



3,(20- iftW[l + i(2S + ^?Si)] 



for Q = c and t For the b-quark we get: 



8 V rri^, cos 2 9w cosa < 



92(H) ^ gi (h)[l + \(2<--^l^)) 

In addition to the above effects one has to consider QCD effects. These effects 
renormalize the quark loops as follows |15|: 



fQCD{q) = h+S.) 

with 

/3(a s ) = ^[l + fa s vr] , 7m (a s ) = 2f 



Thus 



fQCD(q) = l + ^ 



4 7T 

The QCD correction associated with the s-quark loops is: 

f Q CD(q) = l + f t 

The above corrections depend on Q since a s must be evaluated at the scale of mg. 

It convenient to introduce the factor /qcd^) / fQCD(q) into the factors ffi(/i) 
and 92(H) and the factor of /qcd(q') into the the quantities /q. If, however, one 
restricts oneself to the large tanfi regime, the corrections due to the s-quark loops 
is independent of the parameters a and (3 and significant only for the t-quark. 

For large tan (3 we find: 

f c = 0.060 x 1.068 = 0.064, f t = 0.060 x 2.048 = 0.123, f b = 0.060 x 1.174 = 0.070 

quad (model B) 

f c = 0.040 x 1.068 = 0.043, f t = 0.040 x 2.048 = 0.082, f b = 0.040 x 1.174 = 0.047 

quad (model B) 



For a more detailed discussion we refer the reader to Refs. [14, 15 1 
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Table 2: The coupling constants entering £ e //, Eqs. 
solutions #1 - #3. 



and (HD of the text, for 



Quantity 



/3/£ Q 
fvl fv 

f%(H) (model A) 
f l s /f s (model A) 
f s {H) (model B) 
fl(H) (model C) 
fl(NQM) 
f° A (EMC) 

A 



Solution #1 Solution #2 Solution #3 



1.746 x 10 


-5 


2.617 x 10 


-5 


2.864 x 10 


-5 


-0.153 




-0.113 




-0.251 




1.31 x 10~ 


5 


1.30 x 10- 


4 


1.38 x 10- 


-5 


-0.275 




-0.107 




-0.246 




5.29 x 10" 


4 


7.84 x 10- 


3 


6.28 x 10- 


-4 


7.57 x 10~ 


-4 


7.44 x 10- 


3 


7.94 x 10- 


-4 


0.510 x 10 


-2 


3.55 x 10- 


2 


0.704 x 10 


-2 


0.612 x 10 


-3 


0.426 x 10 


-2 


0.844 x 10 


-3 


1.55 x 10~ 


2 


5.31 x 10- 


2 


3.00 x 10- 


-2 



5 Summary of the Input Parameters. 

We have seen that, the vector and axial vector form factors can arise out of Z- 



exchange and s-quark exchange. [12]-|l(| They have uncertainties in them. Here we 



consider the three choices in the allowed parameter space of Kane et al plfl and the 



eight parameter choices of Castano et al [21] These involve universal soft breaking 



masses at the scale. Non-universal masses have also recently been employed [14] 



(see also Arnowitt and Nath [22] In our choice of the parameters the LSP is mostly 
a gaugino. Thus, the Z- contribution is small. It may become dominant in models 
in which the LSP happens to be primarily a Higgsino. Such models, however, are 
excluded by the cosmological bounds on the relic abundance of LSP. The transition 
from the quark to the nucleon level is pretty straightforward in the case of vector 
current contribution. We will see later that, due to the Majorana nature of the 
LSP, the contribution of the vector current, which can lead to a coherent effect 



of all nucleons, is suppressed. [12] The vector current is effectively multiplied by a 
factor of (3 = v/c, v is the velocity of LSP (see Tables 2, 3 ). Thus, the axial current, 
especially in the case of light and medium mass nuclei, cannot be ignored. 

For the isovector axial current one is pretty confident about how to go from the 
quark to the nucleon level. We know from ordinary weak decays that the coupling 
merely gets renormalized from g A = 1 to qa = 1-24. For the isoscalar axial current 
the situation is not completely clear. The naive quark model (NQM) would give 
a renormalization parameter of unity (the same as the isovector vector current). 
This point of view has, however, changed in recent years due to the so-called spin 



crisis, [4C]-[42| i.e. the fact that in the EMC data |4(| it appears that only a 
small fraction of the proton spin arises from the quarks. Thus, one may have to 
renormalize f A by g A = 0.28, for u and d quarks, and g A = —0.16 for the strange 
quarks, jll], |4^[ i.e. a total factor of 0.12. These two possibilities, labeled as NQM 
and EMC, are listed in Tables 2,3. One cannot completely rule out the possibility 



that the actual value maybe anywhere in the above mentioned region. [42] 

The scalar form factors arise out of the Higgs exchange or via s-quark exchange 
when there is mixing [|l4|] between s-quarks qL and qR (the partners of the left- 



12 



Table 3: The same as in Table 2 for solutions #4 - #9. 



Solution 



fvl fv 
f°s(A) 

f°s(B) 

f°s(C) 

fs( A ) 
f° A (NQM) 

f° A (EMC) 
/A 



#4 

0.225 1(T 4 
-0.0809 
-0.179 10~ 4 
-0.531 10~ 2 
-0.315 10~ 2 
-0.207 10~ 5 
6.950 10~ 3 
0.834 10~ 3 
2.490 10~ 2 



_#5 

0.190 10~ 4 
-0.0050 
-0.236 10~ 4 
-0.145 10~ 2 
-0.134 10~ 2 
-0.407 10~ 5 
5.800 10~ 3 
0.696 10~ 3 
1.700 10~ 2 





0.358 10~ 4 
-0.0320 
-0.453 10~ 4 
-0.281 10~ 2 
-0.261 10~ 2 
0.116 10~ 4 
1.220 10~ 2 
0.146 10~ 2 
3.440 10~ 2 



_#7 

0.108 10" 4 
-0.0538 
-0.266 10~ 4 
-0.132 10" 1 
-0.153 10" 1 
0.550 10- 4 
3.760 10~ 2 
0.451 10~ 2 
2.790 10- 1 



#8 

0.694 10~ 4 
-0.0464 
-0.210 10~ 3 
-0.117 10" 1 
-0.118 10" 1 
0.307 10~ 4 
3.410 10~ 2 
0.409 10~ 2 
1.800 10- 1 



J9 

0.864 10~ 4 
-0.0369 
-0.131 10" 3 
-0.490 10~ 2 
-0.159 10~ 2 
0.365 10~ 4 
2.360 10~ 2 
0.283 10~ 2 
2.100 10" 1 



handed and right-handed quarks). We have seen 12] that they have two types of 
uncertainties in them. One, which is the most important, at the quark level due to 
the uncertainties in the Higgs sector. The actual values of the parameters f s and fg 
used here, arising mainly from Higgs exchange, were obtained by considering 1-loop 
corrections in the Higgs sector. As a result, the lightest Higgs mass is now a bit 
higher, i.e. more massive than the value of the Z-boson. [43, Q 



The other type of uncertainty is related to the step going from the quark to 
the nucleon level |Q (see sect. 3.3). Such couplings are proportional to the quark 
masses, and hence sensitive to the small admixtures of qq (q other than u and d) 
present in the nucleon. Again values of f$ and f s in the allowed SUSY parameter 
space are considered (see Tables 2,3). 

6 Expressions for the Unconvoluted Event Rates. 



Combining for results of the previous section we can write 



where 
with 



and 



Jv 



J eff 



Gf 
'y/2 



{(xi7 A 75Xi)<A + (XiXi)^} 



Jx = Nlx{$ + fH + f° Al5 + fh 5 r 3 )N 
= f v {Z) + f v {ql /l = /l(Z) + ft(q) 



f A = f A {Z) + f A {q), f\ = f\(Z) + f\(q) 
J = N{f s +f l s n)N 



(46) 
(47) 

(48) 
(49) 
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We have neglected the uninteresting pseudoscalar and tensor currents. Note 
that, due to the Majorana nature of the LSP, Xi7 A Xi = (identically). 

With the above ingredients the differential cross section can be cast in the form 



2 



diii u 

da(u,v) = — bv)2 [{^S + -g) F\u) + £ flpin Fn(u)] (50) 

ts = {A * m - f s ^n a 4,^r (5D 

C^n = ^^[(4^o(0)) 2 ^ + 2^^o(0)O l( 0)flM + ni(0 ) )a ] (53) 

f A 

Zv = < x o Cv (54) 

_ Qy/mjy) 2 2 /i 2 2 1 2n + l (2u) ] 

^ v ~ fi — V 1 — To" ^ — I \\ — ) 1 ~~ To — rxITT - ! — \2l — 2T v 00 J 

(1 + 4 1 ) 2 /y ^ c (2^ r 6) 2 (1 +7/)^ ( ^ ) 



a 1 o = proton cross-section, i = S, spin, V given by: 

^px = 170 (fs) 2 ( sca l ar ) > (t ne isovector scalar is negligible, i.e. = erf) 

= a ° 3 + /A) 2 (spin) , < x o = a (/° + /* ) 2 (vector) 
where m p is the proton mass, rj = m x jm^A, and p r is the reduced mass and 

(T = —{G F m N f ~ 0.77 x 10~ 38 cm 2 (56) 

Z7T 

u = q 2 b 2 /2 (57) 

or equivalent ly 

« = Q "=A^P (58) 

where b is (the harmonic oscillator) size parameter, q is the momentum transfer to 
the nucleus, and Q is the energy transfer to the nucleus (see Table 4) 
In the above expressions F(u) is the nuclear form factor and 



are the spin form factors |H| (p, p are isospin indices) 
Both form factors are normalized to one at u = 0. 

Slo (Oi) are the static isoscalar (isovector) spin matrix elements (see tables 5 and 
6). 

The non-directional event rate is given by: 

dN p(0) m 

R = Rnon-dir = -JT = "j Ct(m, f ) f (60) 

at m x Amiv 

Where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity and m is the detector 
mass 
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Table 4: The quantity go (forward momentum transfer) in units of fm 1 for three values 
of m x and three typical nuclei. In determining q the value (/3 2 ) 1//2 = 1(C 3 was employed. 







Nucleus 


m x = 30.0 GeV 


m x = 100.0 GeV 


m x = 150.0 GeV 


Ca 


0.174 


0.290 


0.321 


Ge 


0.215 


0.425 


0.494 


Pb 


0.267 


0.685 


0.885 



Table 5: Comparison of the static spin matrix elements for three typical nuclei, Pb 
(present calculation) and 73 Ge, 19 F , 23 Na, 29 Si (see Ref. §|). 



Component 


207 pu 


73 Ge 9/2+ 




23 Na 3/2+ 


29 S*7 ^ 
°n/2 + 


n?(o) 

ni(o)n„(o) 
«3(o) 


0.231 
-0.266 
0.305 


1.005 
-1.078 
1.157 


2.807 
2.707 
2.610 


0.346 
0.406 
0.478 


0.220 
-0.214 
0.208 



The differential non-directional rate can be written as 

dR = dR non _ dir = m da(u,v)\v\ (61) 

m x ArriN 

where da(u,v) was given above. 



The directional differential rate [45] in the direction e is given by : 



dR dir = P(0) m v.eH(v.e) J- da(u, v) (62) 
m x ArriN 2tt 

where H the Heaviside step function. The factor of l/2n is introduced, since the 
differential cross section of the last equation is the same with that entering the 
non-directional rate, i.e. after an integration over the azimuthal angle around the 
nuclear momentum has been performed. In other words, crudely speaking, l/(2ir) is 
the suppression factor we expect in the directional rate compared to the usual one. 
The precise suppression factor depends, of course, on the direction of observation. 



7 Convolution of the Event Rates. 



We have seen that the event rate for LSP-nucleus scattering depends on the relative 
LSP-target velocity. In this section we will examine the consequences of the earth's 
revolution around the sun (the effect of its rotation around its axis is expected to 
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Table 6: Ratio of spin contribution ( 207 Pb/ 73 Ge) at the relevant momentum transfer 
with the kinematical factor 1/(1 + 77) 2 , 77 = m/Am^. 



Solution 


#1 


#2 


#3 


#4 


#5 


#6 


#7 


#8 


#9 


m x (GeV) 


126 


27 


102 


80 


124 


58 


34 


35 


50 


NQM 
EMC 


0.834 
0.645 


0.335 
0.345 


0.589 
0.602 


0.394 
0.499 


0.537 
0.602 


0.365 
0.263 


0.346 
0.341 


0.337 
0.383 


0.417 
0.479 



be negligible) i.e. the modulation effect. In practice this has been accomplished 
by assuming a consistent LSP velocity dispersion, such as a Maxwell distribution 
j2j. More recently other non-isothermal approaches, in the context velocity peaks 
and caustic rings, have been proposed, see e.g Sikivie et al plf . Let us begin with 
isothermal models. In the present paper following the work of Drukier, see Ref. 



where 



[32], we will assume that the velocity distribution is only axially symmetric, i.e. of 
the form 

f(v',\) = N(y esc ,X)(^v )- 3 )[f l (v',X)-f 2 (v',v esc ,X)} (63) 

with 

ffn u K) 2 + (i + A)(K) 2 + K) 2 ) 1 

fi{v , A) = exp[{ J (64) 

^0 

f(> * 1 ^c + A(K) 2 + K) 2 ) 1 
J2{v , v esc , A) = exp[ \ J (65) 

v = \J(2/3){v 2 ) = 220Km/s (66) 

i.e. vq is the velocity of the sun around the center of the galaxy. v esc is the 
escape velocity in the gravitational field of the galaxy, v esc = mbKm/s @. In the 
above expressions A is a parameter, which describes the asymmetry and takes values 
between and 1 and N is a proper normalization constant [[Djj . For y esc — > 00 we 
get the simple expression iV -1 = A + 1 

In a recently proposed non-isothermal model one consider the late in-fall of dark 
matter into our galaxy producing flows of caustic rings. In particular the predictions 
of a self-similar model have been put forward as a possible scenario for dark matter 
density- velocity distribution, see e.g. Sikivie et al [^lj. The implications of such 
theoretical predictions and, in particular, the modulation effect are the subject of 
this section. 

Following Sikivie we will consider 2x N caustic rings, (i,n) , i=(+.-) and n=l,2,...N 
(N=20 in the model of Sikivie et al), each of which contributes to the local density 
a fraction p n of the of the summed density p of each of the i = +,—. It contains 
WIMP like particles with velocity y n = (y nx ,y ny ,y nz ) in units of essentially the 
sun's velocity (vq = 220 Km/ s), with respect to the galactic center. The z-axis is 
chosen in the direction of the disc's rotation, i.e. in the direction of the motion of 
the the sun, the y-axis is perpendicular to the plane of the galaxy and the x-axis 
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Table 7: The velocity parameters y' n = v n /v , y nz = y' nz = v n(jt /v Q , y ny = y' ny = 
Vnz/vo, y nx = y' nx = v nr /v and y n = [(y nz - l) 2 + y 2 ny + ?/ 2 J 1/2 . Also shown are the 
quantities: a n , the caustic rind radii, and p n = 4/[E^°=i4]' ( For the other definitions 
see text ). 



11 


a n (Kpc) 


Vn 


Vnz 






Vn 


Pn 


1 


38.0 


2.818 


0.636 


±2.750 


0.000 


2.773 


0.0120 


2 


19.0 


2.568 


1.159 


±2.295 


0.000 


2.301 


0.0301 


3 


13.0 


2.409 


1.591 


±1.773 


0.000 


1.869 


0.0601 


4 


9.7 


2.273 


2.000 


±1.091 


0.000 


1.480 


0.1895 


5 


7.8 


2.182 


2.000 


0.000 


±0.863 


1.321 


0.2767 


6 


6.5 


2.091 


1.614 


0.000 


±1.341 


1.475 


0.0872 


7 


5.6 


2.023 


1.318 


0.000 


±1.500 


1.533 


0.0571 


8 


4.9 


1.955 


1.136 


0.000 


±1.591 


1.597 


0.0421 


9 


4.4 


1.886 


0.977 


0.000 


±1.614 


1.614 


0.0331 


10 


4.0 


1.818 


0.864 


0.000 


±1.614 


1.619 


0.0300 


11 


3.6 


1.723 


0.773 


0.000 


±1.614 


1.630 


0.0271 


12 


3.3 


1.723 


0.682 


0.000 


±1.591 


1.622 


0.0241 


13 


3.1 


1.619 


0.614 


0.000 


±1.568 


1.615 


0.0211 


14 


2.9 


1.636 


0.545 


0.000 


±1.545 


1.611 


0.0180 


15 


2.7 


1.591 


0.500 


0.000 


±1.500 


1.581 


0.0180 


16 


2.5 


1.545 


0.454 


0.000 


±1.477 


1.575 


0.0165 


17 


2.4 


1.500 


0.409 


0.000 


±1.454 


1.570 


0.0150 


18 


2.2 


1.455 


0.386 


0.000 


±1.409 


1.537 


0.0150 


19 


2.1 


1.432 


0.364 


0.000 


±1.386 


1.525 


0.0135 


20 


2.0 


1.409 


0.341 


0.000 


±1.364 


1.515 


0.0135 



is in the radial direction. We caution the reader that these axes are traditionally 
indicated by astronomers as e^, e r , e z ) respectively. The needed quantities are taken 
from the work of Sikivie et al [^TJ (see Table 7), via the definitions 

Vn = v n/v ,y' nz = V nrj) /VQ = y n z,y' n x = Vnr/Vo = Vnx.Vny = V nz /v = y n y,Pn = 

d n /p,p = J2n=i d n and y n = [(y nz - l) 2 + y 2 ny + y 2 J 1/2 (for each flow +.-). This 
leads to a velocity distribution of the form: 

N 

f(v') = £ 8(v - vo y' n ) (67) 

n=l 

The actual situation, of course, could be a combination of an isothermal contribution 
and late infall of dark matter [jig] . In the present treatment we will consider each 
of these distributions separately. For clarity of presentation we will consider each 
case separately. 

Since the axis of the ecliptic ]13| . lies very close to the y, z plane the velocity of 
the earth around the sun is given by 

ve = + ^i = + v i ( sina x — cosa cos^ y + cosa siwy z ) (68) 

where a is the phase of the earth's orbital motion, a = 2n(t — ti)/TE, where t\ is 
around second of June and Tg = lyear. 
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One can now express the above distribution in the laboratory frame pl| by 
writing v = v + ve 



8 Expressions for the Non-directional Differ- 
ential Event Rate 



The mean value of the non-directional event rate of Eq. (|6l| ) , is given by 

idR\ p(0) m r da(u,v) ^ 3 

\ du / m x ArriN J ' du 

The above expression can be more conveniently written as 
where 

&=f _ 



cffi f \v\ da(u,v) 3 



8.1 No velocity Dispersion-The Case of Caustic Rings 

In the case of caustic rings the last expression takes the form 

<^> = a 2 [t s F {u) + + S spin F spm (n)] (72) 

We remind the reader that p was obtain for each type of flow (+ or -), which explains 



the factor of two. In the Sikivie model 31] we have (2p/p(0) = 1.25, i.e. the whole 



dark matter density lies in the form of caustic rings. In a composite model this can 
only be a fraction of the total density. 

The quantities = S,V,spin are given by Eqs. (|5l|)- fl54|). The quantities 
Fq, Fi, F sp in are obtained from the corresponding form factors via the equations 

F k {u) = F 2 (u)^ k (u) { ^^- , k = 0,l (73) 

F sp in(u) =fii(«)f («) (74) 
The functions ^(u) depend on the model. Introducing the parameter 

5 = — = 0.27, (75) 

in the Sikivie model we find 



3 E Pny 2 n {k ' 1] &(^ ~ u)[(y nz - 1 - - »7 cosaf 

71 = 1 

+ (Vny + ^ cos 7 cosa) 2 + (y nx - - sma) 2 ] 1/2 (76) 
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with 

\/2u r bvo ^ ^ 

Combining the above results the non-directional differential rate takes the form 

/^L\ = R^tT(u)[l-cosaH(u)} (78) 
au p(0) 

In the above expressions R is the rate obtained in the conventional approach [^] by 
neglecting the folding with the LSP velocity and the momentum transfer dependence 
of the differential cross section, i.e. by 

R = ^-^^>[E S + t spm + { ^Lt v ] (79) 

m-y ArriN v <r 



where i = S, V, spin have been defined above, see Eqs fl5l| ) - (|54|). 

The factor T(u) takes care of the u-dependence of the unmodulated differential 
rate. It is defined so that 

duT{u) = 1. (80) 



i.e. it is the relative differential rate. u m j n is determined by the energy cutoff due 
to the performance of the detector, i.e 

Qmin / 01 \ 

Umin = —pz— (81) 

while u m ax is determined the via the relations: 

2 2 

u m ax = min(^f,max(^) , n = 1, 2, N) (82) 

On the other hand H (u) gives the energy transfer dependent modulation amplitude 
(relative to the unmodulated amplitude). The quantity t takes care of the modi- 
fication of the total rate due to the nuclear form factor and the folding with the 
LSP velocity distribution. Since the functions Fq(u), F\ and F sp i n have a different 
dependence on u, the functions T(u)and H(u) and t, in principle, depend some- 
what on the SUSY parameters. If, however, we ignore the small vector contribution 
and assume (i) the scalar and axial (spin) dependence on u is the same, as seems 
to be the case for light systems pQ|'p9|, or (ii) only one mechanism (S, V, spin) 
dominates, the parameter R contains the dependence on all SUSY parameters. The 
parameters t and T(u) depend on the LSP mass and the nuclear parameters, while 
the H(u) depends only on the parameter a. 

8.2 Velocity Dispersion-Isothermal Models 

Expanding in powers of 5 , see Eq. (75) and keeping terms up to linear in it we can 
manage to perform the angular integrations |l8| in Eq. (|n]) and get Eq. (|72|). Now 
the quantities ^k( u ) are given by 

= $(p) t ki a Vu) +0.135 cos a^ (1)ifc (a v / n)] (83) 
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and 

$(l),k( x ) = N (Vesc A)e- A (e^ 1 $ (0jfc (x) - exp[-yl sc ]$ m {x)) (84) 
* (0)fc (x) = -|= d^-^pC-Cl + A)y 2 ))(^(A, (A + 1)2!/) + G,(A, y))) (85) 

VD7r 

= jT" dyy 2fe - 1 exp(-Ay 2 ))G;(A,y)) (86) 

In the above expressions 

G (0,y) = , Gi(0,2/)=0 (87) 

F (A, x) = (A + 1)~ 2 x sinh(x) (88) 

Fx (A, a?) = (1 + A)" 2 [(2 + A)(x/2) + 1) [x cosh(x) - (2A + 3) sinh(x)] (89) 

note that here x = (A + l)2y. The functions G cannot be obtained analytically, but 
they can easily be expressed as a rapidly convergent series in y = ^ , which will not 
be given here. 

The functions Gj(A,y), associated with the small second term of the velocity 
distribution are obtained similarly [18]. 



The non-directional differential rate takes the form 



<^> = 5^tr(u)[(l + coBaff(tt))] (90) 



With R given by Eq. (|79| ) and p (0) is the part of the total LSP density attributed to 
this mode. Note the difference of sign in the definition of the modulation amplitude 
H compared to Eq. (|78|) . 

Here u m i n is determined by the energy cutoff due to the performance of the 
detector. u max is determined by the escape velocity v esc via the relation: 

""max — o v 31 -) 

a z 

Considering only the scalar interaction we get R — > Rs and 

tT(u)=a 2 F\u)i, {0)fi (aV^) (92) 

For the spin interaction we get a similar expression except that R — > -R sp j n and 
.F 2 — ► i^ii. Finally for completeness we will consider the less important vector 
contribution. We get R — > -Ry and 

1 2?y + 1 2a 2 
t T{u) = F' ! (n)[^ (0 ) i i(a v / «) - ^ ^ 2 ^ - ^ «^o),o(^]y (93) 

The quantity T(u) depends on the nucleus through the nuclear form factor or the 
spin response function and the parameter a. The modulation amplitude takes the 
form 

ibr-i\ u(a\/u) , , 

H(u) = 0.135 - ( ] ' (94) 
■0(o),fe (aVu) 

Thus in this case the H(u) depends only on a^/u, which coincides with the parameter 



x of Ref. [35], i.e. only on the momentum transfer, the reduced mass and the size 
of the nucleus. 

Returning to the differential rate it is sometimes convenient to use the quantity 
T(u)H(u) rather than H, since H(u) may appear artificially increasing function of 
u due to the faster decrease of T(u) (H(u) was obtained after division by T(u)) 
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9 Expressions for the Directional Differential 
Event Rate 



There are now experiments under way aiming at measuring directional rates [38] 
using TPC counters which permit the observation of the recoiling nucleus is observed 
in a certain direction. Prom a theoretical point of view the directional rates have 



been previously discussed by Spegel |57]] and Copi et al [45|. The rate will depend 
on the direction of observation, showing a strong correlation with the direction of 
the sun's motion. In a favorable situation the rate will merely be suppressed by a 
factor of 2ir relative to the non-directional rate. This is due to the fact that one 
does not now integrate over the azimuthal angle of the nuclear recoiling momentum. 
The directional rate will also show modulation due to the Earth's motion. We will 
again examine an non-isothermal non symmetric case ( the Sikivie model) and a 
3-dimensional Gaussian distribution with only axial symmetry. 



The mean value of the directional differential event rate of Eq. (62), is defined 

by 

/dR\ _ p(0) m 1 r ^ cH(v c) ^ a ( Ul v ^ d z v (95) 
\ du I dir m x Am n 27r J ' du 

where e is the unit vector in the direction of observation. It can be more conveniently 
expressed as 

(f) = PM ^-xf^)(f) 06) 
\du/dir m x Am N v x du dir 

where 

,dH x If v.eH(v .e) e/ ^da(u,v) 3 



< ^ > *=w^r (, '' , '* ) ^ <i " (97) 

9.1 Directional Differential Event Rate in the Case of 
Caustic Rings. 



The model of Sikivie et al |31j| , which is a non isothermal and asymmetric one, 
offers itself as a perfect example for the study of directional rates. So we are going 
to begin our discussion with such a case. Working as in the previous section we get 

@ 

dT, 2d a 2 (v 2 ) - 

{-^)dir = —[t s F (u) + ^Lt v F 1 (u) + t spin F spin (u)] (98) 

where the Sj, i = S, V, spin are given by Eqs. (|5l|)- fl54|) . The quantities Fq, Fx, F spin 
are obtained from the equations 

F k (u) = F 2 (n)^(n)§-^, k = 0, 1 (99) 
F sp in(u) = F n (u)^ (u) (100) 

In the Sikivie model we find 

[2 N v 2 5 

^k{u) = J - Pnyl {k ~ 1] &(^ ~ u)\{y nz - 1 - - sin-/ cosa)e z .e 

V o ^ (X Zi 

n=l 

s & 

+ iVny + x cos l cosa)e y .e + (y nx - - sina)e x .e\ (101) 
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In the model considered here the z-component of the LSP's velocity with respect 
to the galactic center for some rings is smaller than sun's velocity, while for some 
others it is larger. The components in the y and the x directions are opposite for 
the + and - flows. So we will distinguish the following cases: a) e has a component 
in the sun's direction of motion, i.e. < 9 < tt/2, labeled by u (up), b) Detection 
in the direction specified by tt/2 < 9 < tt, labeled by d (down). The differential 
directional rate takes a different form depending on which quadrant the observation 
is made. Thus : 

1. In the first quadrant (azimuthal angle < <j> < tt/2). 

( ^7 } = R W)^ T{u)m{u) ~ c ° s a Hl{u))lez - eI 

+ (Ry + cosaH l 2 (u) H ^—?-(\cosa\ + cosa))\e y .e\ 

+ (Rl - sinaHi(u) + ^^-(\sina\ - sina))\e x .e\] (102) 

2. In the second quadrant (azimuthal angle tt/2 < <fi < tt) 

<^> = R^^T(u)[(Ri(u)-cosaHi(u))\e z .e\ 

+ (R y + cosaH 2 (u) H ^-^-(\cosa\ - cosa))\e y .e\ 

+ (Ri + sinaHi(u) H ^-(\sina\ + sina))\e x .e\] (103) 

3. In the third quadrant (azimuthal angle tt < 4> < 37r/2). 

W = R^^ T (u)[(Rl(u)-co S aHl(u))\e z .e\ 

+ (R y — cosaH l 2 (u) H ^ — (|cosa| — cosa))\e y .e\ 

+ (Ri + sinaHl(u) + ^ ' (\sina\ + sina))\e x .e\] (104) 

4. In the fourth quadrant (azimuthal angle 37r/2 < <p < 2tt) 

<^> = R^^T(u)[(Ri(u)- cos a H{(u))\e z .e\ 

H^(u) 

+ (R y — cosaH' 2 (u) H ^ — (|cosa| — cosa))\e y .e\ 

+ (Ri - sinaHi(u) H s -^-(\sina\ - sma))|e x .e|] (105) 

where i = u,d 

By the reasoning given above, if one mechanism is dominant, the parameters 
R x , R y , R z , Hi, i?2 5 H%, H c , H s for both directions u and d depend only on fi r and a. 
They are all independent of the other SUSY parameters. 
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9.2 The Directional Differential Event Rate in the Case 
of Velocity Dispersion 

The dependence of the rate depends on the direction of observation in a rather 
complicated way. The integrals can only be done numerically [45]. To simplify 
matters we made a power expansion in 6 and kept terms up to linear in it. To make 
the presentation tractable we will will give expressions valid only for directions of 
greatest interest, i.e. close to the coordinate axes. In the sun's direction of motion 
we have a modulated as well as a non modulated amplitude. In the other two 
directions we only have a modulated amplitude. Unlike the case of caustic rings, 
now the direction opposite to the sun's direction of motion is favored. We found 
it more convenient, however, to present our results as the absolute value of the 
difference of the rates in the directions e and — e. 
Working as in the previous subsection we get [18] 

dH 1 - (v 2 ) - 

(^)dir = ^ a 2 [t s F (u) + ^-J-tvF^u) + E spin F spin (u)} (106) 
The quantities Fo , F\ , F sp i n are obtained from the equations 

F k {u) = F 2 (u)^ k (u)^^f, k = 0,l (107) 



Now 



1 



^k{u) = 2[(ip(o),k( a ^ + 0.135 cos a^(i),fc(aVM))| 



e z .e 



— 0.117cos aV'(2), k( a Vu)\ e y e \ + 0.135 sin a^(^\ k (a-y/u)\e x .e\] (108) 

with 

i) m {x) = N(y esc , A)e" A (e- 1 $ (0>fc (x) - exp[-y 2 esc }<S>' m (x)) (109) 

rVeSC 2k-l < ( 1 , 



* (0>fc (x) = -= / dyy 2k -'exp(-(l + X)y 2 ))(F l (X,2(X + l)y) + G l (X,y))) (110) 



*' w>Jfc (s) = -|= J^dyy^expi-Xy^G'^y)) (111) 
In the above expressions 

Fi{X, x) = xcoshx - sinhx , i = 0,2,3 (112) 

Fi(A,x) = 2(1 -A) [( %^%- +l)sinh X - xcoshx] (113) 

The purely asymmetric quantities Gi satisfy 

Gi(0,y) = 0, » = 0,4 (114) 

the qualities G^(0, y) = 0, i = 0, 4 refer to the second term of the Velocity distribu- 
tion and were obtained in an analogous fashion. 

If we consider each mode (scalar, spin vector) separately the directional rate 
takes the form 



dR - p(0) t° R° 

d^ Uir ~ p(0) 4vr 

(115) 



(-j-)dir = R — — | (1 + cos aHi(u))e z .e — cosaH 2 (u)e y .e + sinaH^(u)e x .e 
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In other words the directional differential modulated amplitude is described in 
terms of the three parameters, Hi(u), 1=1,2 and 3. The unmodulated amplitude 
R°(u) is again normalized to unity. The parameter t° entering Eq. (|115| ) takes care 
of whatever modifications are needed due to the convolution with the LSP velocity 
distribution in the presence of the nuclear form factors. 



From Eqs. (|106j) - (115) we see that if we consider each mode separately the 



differential modulation amplitudes H (I) take the form 

*,(.)- 0.135^|^ . 1-1,3 , J7 >M _ail73$£sg (116) 

Thus in this case the Hi depend only on ay/u, which coincides with the parameter x 
of Ref. [ j33|| . This means that Hi essentially depend only on the momentum transfer, 
the reduced mass and the size of the nucleus. We note that in the case A = we 
have H 2 = 0.117 and H 3 = 0.135 

It is sometimes convenient to use the quantity Ri rather than Hi defined by 

Ri = R H h I = 1,2,3. (117) 

The reason is that Hi, being the ratio of two quantities, may appear superficially 
large due to the denominator becoming small. 



10 The Total Non-directional Event Rates 

Integrating Eq. (|7^) we obtain for the total non-directional rate in the case of 
caustic rings the expression: 

R = R -77^ t [1 - h(a, Q min )cosa)] (118) 

where Qmin is the energy transfer cutoff imposed by the detector. The modulation 
is described by the parameter h. Similarly integrating Eq. (|x]) we obtain for the 
total non-directional rate in our isothermal model as follows: 

R = R^t[(l + h{a,Q min )cosa)} (119) 

Note the difference of sign in the definition of the modulation amplitude h compared 
to Eq. ( PSD , where Q min is the energy transfer cutoff imposed by the detector. 
The modulation can be described in terms of the parameter h. 

The effect of folding with LSP velocity on the total rate is taken into account via 
the quantity t. The SUSY parameters have been absorbed in R. From our discussion 
in the case of differential rate it is clear that strictly speaking the quantities t and 
h also depend on the SUSY parameters. They do not depend on them, however, if 
one considers the scalar, spin etc. modes separately. 

The meaning of t is clear from the above discussion. We only like to stress 
that it is a common practice to extract the LSP nucleon cross section from the the 
expected experimental event rates in order to compare with the SUSY predictions 
as a function of the LSP mass. In such analysis the factor t is omitted. It is clear, 
however, that, in going from the data to the cross section, one should divide by t. 
The results will be greatly affected for large reduced mass. 
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11 The Total Directional Event Rates 



We will again examine separately the case of caustic rings and the isothermal models 
considered above. 



11.1 The Total Directional Event Rates in the Case of 
Caustic Rings 



Integrating Eqs. ( |102| ) - ( |105| ) we obtain: 

1. In the first quadrant (azimuthal angle < (ft < n/2) 

Rdir = Rj^T^KK-cosa h\)\e z .e\ 



+ (r* + cosah l 2 + -^-(\cosa\ + cosa))\e y .e\ 

hi 



+ (r x — sinah\ + -^-(\sina\ — sina))\e x .e\] (120) 



2. In the second quadrant (azimuthal angle ir/2 < (ft < it) 

R * ir = ^o)i [(r * ~ cosa h ^ Gz - e \ 

hi 



+ (r l y + cosah 2 (u) + -^-(\cosa\ — cosa))\e y .e\ 
h* 

+ (r x + sinahl + -^-(\sina\ + sina))\e x .e\] (121) 

3. In the third quadrant (azimuthal angle ir < (ft < 3tt/2). 

Rdir = Rj^^[(ri-cosa h\)\e z .e\ 

h} (u) 

+ [r l y — cosah l 2 (u) H — 2- — (|cosa| — cosa))\e y .e| 
h} 

+ (r* +sinaH\ + -^-(\sina\ + sina))\e x .e\] (122) 

4. In the fourth quadrant (azimuthal angle 3ir/2 < (ft < 2ir) 

Rdir = R^7^[(<-cosa h\)\e z .e\ 

K 

+ (r y — cosah 2 + -^-(\cosa\ — cosa))\e y .e\ 
W 

+ (r x — sinah\ + -^-{\sina\ — sina))\e x .e\] (123) 

11.2 The Total Directional Event Rates in Isothermal 
Models 

We remind the reader that in this case we take the difference of the rates in two 
opposite directions. 
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Integrating Eq. (99) we obtain 



Rdir = R[p\0)/(p(0))(t°/4Tr)\(l + h 1 {a,Q min )cosa)e z .e 

- h 2 (a, Qmin) cosae y .e + h 3 (a, Q m in) sinae x .e\ (124) 

note that in the above expressions, unlike the case of caustic rings, the rate is 
normalized to t /2 and not to t In other words the effect of folding with LSP 
velocity on the total rate is taken into account via the quantity t°. All other SUSY 
parameters have been absorbed in R, under the assumptions discussed above in the 
case of non-directional rates. 

We see that the modulation of the directional total event rate can be described 
in terms of three parameters hi, 1=1,2,3. In the special case of A = we essentially 
have one parameter, namely hi, since then we have h 2 = 0.117 and /13 = 0.135. 

Given the functions hi(a, Qmin) one can plot the the expression in Eq. ( |124| ) as 
a function of the phase of the earth a. 



12 Results and Discussion 



The three basic ingredients of our calculation were the input SUSY parameters (see 
sect. 5), a quark model for the nucleon (see sect. 4) and the structure of the nuclei 
involved (see sect. 6). The input SUSY parameters used for the results presented in 
Tables 1,2 and 3 have been calculated in a phenomenologically allowed parameter 
space (cases #1, #2, #3 of Kane et al [^] and cases #4-9 of Castano et al [21]. 
Our own SUSY parameters will appear elsewhere [19]. 

For the coherent part (scalar and vector) we used realistic nuclear form factors 
and studied three nuclei, representatives of the light, medium and heavy nuclear 
isotopes (Ca, Ge and Pb). In Tables 8,9 and 10 we show the results obtained 
for three different quark models denoted by A (only quarks u and d) and B, C 
(heavy quarks in the nucleon). We see that the results vary substantially and 
are very sensitive to the presence of quarks other than u and d into the nucleon. 
The spin contribution, arising from the axial current, was computed in the case 
of a number of both light and heavy nuclei, including the 207 Pb system. For the 
isovector axial coupling the transition from the quark to the nucleon level is trivial (a 
factor of qa = 1-25). For the isoscalar axial current we considered two possibilities 
depending on the portion of the nucleon spin which is attributed to the quarks, 
indicated by EMC and NQM. [13] The ground state wave function of 20S Pb was 
obtained by diagonalizing the nuclear Hamiltonian f4^[ - ||48| in a 2h-lp space which 
is standard for this doubly magic nucleus. The momentum dependence of the matrix 
elements was taken into account and all relevant multipoles were retained (here only 
monopole and quadrupole). 

In Table 5, we compare the spin matrix elements at q = for the most popular 
targets considered for LSP detection 207 Pb, 73 Ge, 19 F 23 Na and 29 Si. We see 
that, even though the spin matrix elements Q 2 in the case of 207 Pb are about 
factor of three smaller than those for 73 'Ge obtained in Ref. [^] (see Table 5), their 
contribution to the total cross section is almost the same (see Table 6) for LSP 
masses around 100 GeV. Our final results for the quark models (A, B, C, NQM, 
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Table 8: The quantity {dN/dt)o = Rt in y~ l Kg~ 1 and the modulation parameter h for 
the coherent vector and scalar contributions in the cases #1 - #3 and for three typical 
nuclei. 





Vector Contribution 


Scalar Contribution 




{dN/dt) h 


(dN/dt} h 


Case 


(xicr 3 ) 


Model A Model B Model C 


#1 
Pb #2 
#3 


0.264 0.029 
0.162 0.039 
0.895 0.038 


0.151 x 10~ 3 0.220 0.450 -0.002 
0.410 x 10" 1 142.860 128.660 0.026 
0.200 x 10~ 3 0.377 0.602 -0.001 


#1 

Ge #2 
#3 


0.151 0.043 
0.053 0.057 
0.481 0.045 


0.779 x 10" 4 0.120 0.245 0.017 
0.146 x 10" 1 51.724 46.580 0.041 
0.101 x 10" 3 0.198 0.316 0.020 


#1 

Ca #2 

#3 


0.079 0.053 
0.264 0.060 
0.241 0.053 


0.340 x 10~ 4 0.055 0.114 0.037 
0.612 x 10~ 2 22.271 20.056 0.048 
0.435 x 10~ 4 0.090 0.144 0.038 



EMC) are presented in Tables 8, 9 for SUSY models #l-#3 of Kane et al [21] and 
Table 10 for SUSY models #4-#9 of Castano et al ^ 

In discussing the effects of folding with the LSP velocity combined with the 
nuclear form factor and specialized our results for the target 127 1. To this end we 
considered only the scalar interaction and studied the effects of the detector energy 
cutoffs, by considering two typical cases Q m %n = 10, 20 KeV. 

Special attention was paid to the the directional rates and the modulation effect 
due to the annual motion of the earth. 

We will start our discussion with the non isothermal velocity spectrum due to 
caustic rings resulting from the self- similar model of Sikivie et al [31]. 

The total rates are described in terms of the quantities t, r l x , r l y , r\ for the unmod- 
ulated amplitude and h,h\,h % 2,h l s ,h l c ,h l s i = u,d for the modulated one. In Table 
11 I we show how these quantities vary with the detector energy cutoff and the LSP 
mass. Of the above list only the quantities t and h enter the non-directional rate. 
We notice that the usual modulation amplitude h is smaller than the one arising 
in isothermal models |l7], 18[. The reason is that there are cancelations among the 
various rings, since some rings are characterized by y nz > 1, while for some others 
Vnz < 1 (see Table 7). Such cancelations are less pronounced in the isothermal 
models. As expected, the parameter t, which contains the effect of the nuclear form 
factor and the LSP velocity dependence, decreases as the reduced mass increases. 

In the case of isothermal models we will limit ourselves to the discussion of the 
directional rates. In the special case of the direction of observation being close to 
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Table 9: The spin contribution in the LSP — 207 Pb scattering for two cases: EMC data 
and NQM Model for solutions #1, #2, #3. 





EMC DATA 


NQM MODEL 


Solution 


(dN/dt) {y- l Kg- 1 ) 


h 


(dN/dt) {y- l Kg- 1 ) h 


#1 


0.285 x 10 -2 


0.014 


0.137 x 10~ 2 0.015 


#2 


0.041 


0.046 


0.384 x 10~ 2 0.056 


#3 


0.012 


0.016 


0.764 x 10" 2 0.017 



Table 10: The same quantities as in Table 8 in the case of Pb for the solutions #4 — #9. 
#8 snd #9 are no-scale models. The values of (dN/dt) = Rt for Model A and the Vector 
part must be multiplied by xlCT 2 . 



Scalar Part 



Vector Part 



Spin Part 



(dNx 
\ dt /o 



ldN\ 
\ dt /0 



/ dN\ 
j dt /0 



Case 



B 



C 



EMC 



NQM 



#4 
#5 
#6 
#7 
#8 
#9 



0.03 
0.46 
0.16 
4.30 
2.90 
2.90 



22.9 
1.8 
5.7 
110.0 

73.1 
1.6 



8.5 
1.4 
4.8 
135.0 
79.8 
1.7 



0.003 
-0.003 
0.007 
0.020 
0.020 
0.009 



0.04 
0.03 
0.11 
0.94 
0.40 
0.95 



0.054 
0.053 
0.057 
0.065 
0.065 
0.059 



0.80 10" 3 
0.37 10" 3 
0.44 10" 3 
0.67 
0.22 
0.29 



0.16 10" 2 
0.91 10" 3 
0.11 10" 2 
0.87 
0.35 
0.37 



0.015 
0.014 
0.033 
0.055 
0.055 
0.035 



the coordinate axes the rate is described in terms of the three quantities to an d 
hi, i = 1,2,3 (see Eq. ( |124| )). These are shown in tables 12-14 for various values 
of Qmin and A. For the differential rate the reader is referred to our previous work 
p7| , |lg| ]. We mention again that /12 and h% are constant, 0.117 and 0.135 respectively, 
in the symmetric case. On the other hand hi,h,2 and /13 substantially increase in 
the presence of asymmetry. 



13 Conclusions 



In the present paper we have calculated the parameters, which described the event 
rates for direct detection of SUSY dark matter. We found that the event rates are 
quite small and only in a small segment of the allowed parameter space they are 
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Table 11: The quantities t and h entering the total non-directional rate in the case of 
the target 53/ 127 for various LSP masses and Q m in in KeV. Also shown are the quantities 
r l j, hj i — u,d and j = x, y, z, c, s, entering the directional rate for no energy cutoff. For 
definitions see text. 











t 

IjO-T 


mass 


in PoV 
111 uc V 






Quantity 


Martin 


1 o 

±U 


in 

oU 


OU 


oU 


1 on 

1UU 


1 OK 
1ZO 


9^n 
zOU 


t 


n 
u.u 


1.401 


1 079 
1 .u I z 


O 7^1 
U. / 01 


n a 11 

U.4 ( ( 


n 170 
u.o / y 


O QOQ 
U.0U0 


O 1 71 

u.i 1 


In 

n 


n n 
u.u 


099 

u.uzz 


091 

U.UZO 


O 09/1 
U.Uz4 


n 09^. 

U.UZD 


O 09K 
U.UzO 


O 09fi 
U.UZO 


09fi 

u.uzo 


r z 


n 
u.u 


o 79fi 
U. ( zo 


717 

U. ( ( 


O 7/1 7 
U. f 4 f 


o iki 

U. f O f 


n 7ko 
u. / OU 


O 7fi1 
U. i Dl 


O 7K1 
U. 1 01 


r y 


n n 
u.u 


o 9/1^ 

U.Z40 


U.ZOl 


O 91 Q 

u.ziy 


911 
U.Z1 ± 


O 90Q 

u.zuy 


O 90S 
U.ZUo 


908 
U.ZUo 


r x 


n 
u.u 


o 


n 1^1 

U.oOl 


O IRK 
U.oDD 


*377 


O ISO 
U.ooU 


O ^&1 


181 
U.ool 


h u 

n z 


n 
u.u 


o 09£ 
u.uzo 


097 

u.uz < 


O 098 
U.Uzo 


090 
U.UZy 


O 09Q 

u.uzy 


O OQO 
U.UoU 


OIO 
U.UoU 


n y 


n 
u.u 


n 09i 

U.UZ1 


n 091 

U.UZ1 


090 

u.uzu 


090 
u.uzu 


01 Q 

u.uiy 


01 Q 

u.uiy 


01 Q 

u.uiy 


uu 


n n 


041 


0.044 


046 


U.U4:U 


048 


04Q 


049 


K 


0.0 


0.036 


0.038 


0.040 


0.041 


0.042 


0.042 


0.042 


uu 
n s 


n 

u.u 


n oifi 

u.uou 


094 
U.UZ^t 


094 


091 


091 


099 
u. uzz 


099 
u.uzz 


r d 
1 z 


0.0 


0.274 


0.263 


0.253 


0.243 


0.240 


0.239 


0.239 


r d 

1 y 


0.0 


0.019 


0.011 


0.008 


0.007 


0.007 


0.007 


0.007 


r x 


o n 
u.u 


n 9/L^ 
u.z^o 


941 
U.Z40 


9^fi 


997 
U.ZZ ( 


99^ 
u.zzo 


99^ 
u.zzo 


99^ 
U.ZZO 


hi 


0.0 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


0.004 


K 


0.0 


0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


K 


0.0 


0.022 


0.021 


0.021 


0.020 


0.020 


0.020 


0.020 


h d c 


0.0 


0.019 


0.018 


0.018 


0.017 


0.017 


0.017 


0.017 


hi 


0.0 


0.001 


0.001 


0.000 


0.000 


0.000 


0.000 


0.000 


t 


10.0 


0.000 


0.226 


0.356 


0.265 


0.224 


0.172 


0.098 


h 


10.0 


0.000 


0.013 


0.023 


0.025 


0.025 


0.026 


0.026 


t 


20.0 


0.000 


0.013 


0.126 


0.139 


0.116 


0.095 


0.054 


h 


20.0 


0.000 


0.005 


0.017 


0.024 


0.025 


0.026 


0.026 



above the present experimental goals. One, therefore, is looking for characteristic 
signatures, which will aid the experimentalists in reducing background. These are 
two: a) The non directional event rates, which are correlated with the motion of the 
Earth (modulation effect) and b) The directional event rates, which are correlated 
with both the velocity of the sun, around the center of the Galaxy and the velocity 
of the Earth. We separated our discussion into two parts. The first deals with the 
elementary aspects, SUSY parameters and nucleon structure, and is given in terms 
of the nucleon cross-section. The second deals with the transition from the nucleon 
to the nuclear level. In the second step we also studied the dependence of the rates 
on the energy cut off imposed by the detector. We presented our results in a fashion 
understandable by the experimentalists. We specialized our results in the case of 
the coherent process in the case of 127 1, but we expect the conclusions to be quite 
general. 

The needed local density and velocity spectrum of the LSP were obtained in 
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Table 12: The quantities t°, hi and h m for A = in the case of the target 53/ for various 
LSP masses and Q m in in KeV (for definitions see text). Only the scalar contribution is 
considered. Note that in this case h 2 and h% are constants equal to 0.117 and 0.135 
respectively. 











LSP 


mass 


in GeV 






Quantity 


Qmin 


10 


30 


50 


80 


100 


125 


250 


t° 


0.0 


1.960 


1.355 


0.886 


0.552 


0.442 


0.360 


0.212 


hi 


0.0 


0.059 


0.048 


0.037 


0.029 


0.027 


0.025 


0.023 


t° 


10. 


0.000 


0.365 


0.383 


0.280 


0.233 


0.194 


0.119 


hi 


10. 


0.000 


0.086 


0.054 


0.038 


0.033 


0.030 


0.025 


t° 


20. 


0.000 


0.080 


0.153 


0.136 


0.11 


0.102 


0.065 


hi 


20. 


0.000 


0.123 


0.073 


0.048 


0.041 


0.036 


0.028 



two special classes: 1) Non isothermal models and 2) Isothermal models. As we 
have already mentioned the actual situation may be a combination of an isothermal 
contribution and late infall of dark matter. In the present treatment we consider 
each of these distributions separately. 

In the first case case we assumed a late in-fall of dark matter into our galaxy, 
the needed parameters were taken from the work of Sikivie et al [31] in the context 
of a self-similar model, which yields 40 caustic rings. Our results, in particular the 
parameters t, see Table 11, indicate that for large reduced mass, the kinematical 
advantage of fi r (see Eqs. (pi]) - (|54|) is partly lost when the nuclear form factor 
and the convolution with the velocity distribution are taken into account. Also, if 
one attempts to extract the LSP-nucleon cross section from the data, in order to 
compare with the predictions of SUSY models, one must take t into account, since 
for large reduced mass t is different from unity. 

In the case of the non-directional total event rates we find that the maximum 
no longer occurs around June 2nd, but about six months later. The difference 
between the maximum and the minimum is about 4%, smaller than that predicted 
by the symmetric isothermal models [17, 18|. In the case of the directional rate we 
found that the rates depend on the direction of observation. The biggest rates are 
obtained, If the observation is made close to the direction of the sun's motion. The 
directional rates are suppressed compared to the usual non-directional rates by the 
factor f re( i = n/(2ir). We find that k = r" ~ 0.7, if the observation is made in the 
sun's direction of motion, while k ~ 0.3 in the opposite direction. The modulation 
is a bit larger than in the non-directional case, but the largest value, 8%, is not 
obtained along the sun's direction of motion, but in the x-direction (galactocentric 
direction). 

In the case of the isothermal models we restricted our discussion to the direc- 
tional event rates. The reduction factor of the total directional rate, along the sun's 
direction of motion, compared to the total non directional rate depends, of course, 
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Table 13: The same as in the previous table, but for the value of the asymmetry parameter 
A = 0.5. 











LSP 


mass 


in GeV 






Quantity 


O 

Trim 


10 


30 


50 


80 


100 


125 


250 


t° 


0.0 


2.309 


1.682 


1.153 


0.737 


0.595 


0.485 


0.288 


hi 


0.0 


0.138 


0.128 


0.117 


0.108 


0.105 


0.103 


0.100 


h 2 


0.0 


0.139 


0.137 


0.135 


0.133 


0.133 


0.133 


0.132 


h* 


0.0 


0.175 


0.171 


0.167 


0.165 


0.163 


0.162 


0.162 


t° 


10. 


0.000 


0.376 


0.468 


0.365 


0.308 


0.259 


0.160 


hi 


10. 


0.000 


0.174 


0.139 


0.120 


0.114 


0.110 


0.103 


h 2 


10. 


0.000 


0.145 


0.138 


0.135 


0.134 


0.134 


0.133 


h 3 


10. 


0.000 


0.188 


0.174 


0.167 


0.165 


0.164 


0.162 


t° 


20. 


0.000 


0.063 


0.170 


0.171 


0.153 


0.134 


0.087 


hi 


20. 


0.000 


0.216 


0.162 


0.133 


0.124 


0.118 


0.107 


h 2 


20. 


0.000 


0.155 


0.143 


0.137 


0.136 


0.135 


0.133 


h 3 


20. 


0.000 


0.209 


0.182 


0.171 


0.168 


0.166 


0.164 



on the nuclear parameters, the reduced mass and the asymmetry parameter A [18]. 



It is given by the parameter f re d = to/(4n t) = k/(2it). We find that k is around 
0.6 for no asymmetry and around 0.7 for maximum asymmetry (A = 1.0). In other 
words it is not very different from the naively expected f re d = l/(27r), i.e. k = 1. 
The modulation of the directional rate increases with the asymmetry parameter A 
and it depends, of course, on the direction of observation. For Q m %n = it can reach 
values up to 23%. Values up to 35% are possible for large values of Q m0 n, but they 
occur at the expense of the total number of counts. In all cases our results, in partic- 
ular the parameters t, see Table 11, and to, see Tables 12-14, indicate that for large 
reduced mass, the kinematical advantage of fj, r (see Eqs. (51)- (|3) is partly lost 



when the nuclear form factor and the convolution with the velocity distribution are 
taken into account. To be more precise, if one attempts to extract the LSP-nucleon 
cross section from the data, in order to compare with the predictions of SUSY mod- 
els, one must take t into account, since for large reduced mass t is different from 
unity. Acknowledgments: Partial support by TMR No ERB FMAX-CT96-0090 of 

the European Union is happily acknowledged. 
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